Persistence and eventual demise of oxygen molecules at terapascal pressures 
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Computational searches for structures of solid oxygen under pressures in the multi TPa range 
have been carried out using density-functional-theory methods. We find that molecular oxygen 
persists to about 1.9 TPa at which it transforms into a semiconducting square spiral-like polymeric 
structure (/4i /acd) with a band gap of ~3.0 eV. Solid oxygen forms a metallic zig-zag chain-like 
structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure 
is increased and a structure of Fmmm symmetry forms at about 9.5 TPa in which each atom has 
four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with 
compression, although the structure becomes more symmetric. The electronic properties of oxygen 
have a complex evolution with pressure, swapping between insulating, semiconducting and metallic. 



PACS numbers: 61.50.Ks 71.20.-b 81.05.Zx 81.40.Vw 

Among the simple molecules studied at high pressures, 
oxygen has attracted particular attention due to its fun- 
damental importance and intriguing properties. ^{5] For 
example, oxygen is the only known elemental molecule 
that exhibits magnetism, which adds substantial com- 
plexity to its phase diagram. When liquid oxygen is 
cooled at ambient pressure it undergoes a sequence of 
transitions to the 7, (3 and a solid phases at 54.39, 43.76 
and 23.88 K, respectively [6]. Upon increase of pres- 
sure the monoclinic a-02 (C2/m) phase transforms into 
an orthorhombic S-O2 {Fmmm) phase at about 3 GPa 
and to e-02 at about 10 GPa. The structure of e-02 
has only recently been solved by x-ray diffraction (XRD) 
studies of single crystal [3] and powder [4] samples, al- 
though its vibrational spectrum was reported more than 
20 years earlier [7]. The unit cell of e-02 has C2/m sym- 
metry and contains four O2 molecules forming Og units. 
Both a and 6-O2 are antiferromagnetic and the mag- 
netic collapse at the ^-e transition was predicted using 
molecular dynamics simulations [8]. This breakdown of 
the long-range antiferromagnetic order at about 8 GPa 
was recently observed in a neutron scattering experiment 
[2]. Density functional theory (DFT) studies have found 
another chain-like structure to be energetically slightly 
more favorable than the Os structure [9l [10] , although it 
has not been observed in experiments. 

Recent experiments [5] have shown that insulating e-02 
remains stable up to about 96 GPa before undergoing a 
continuous displacive and isosymmetric transition to the 
C phase, in agreement with earlier predictions [10] . The 
metallic ^ phase [11] [H] has C2/m symmetry [2] and 
superconducts at temperatures below 0.6 K [T]. Many- 
body perturbation theory GW calculations [13l [M] have 
suggested that the metal-insulator transition occurs at 
lower pressures than the measured e-C transition pressure 



of 96 GPa. The above information, however, covers only 
a small part of the phase diagram and rather little is 
known about pure oxygen at pressures above 100 GPa. 

A previous DFT study reported that oxygen molecules 
persist to at least 250 GPa [10]. It is interesting to 
speculate about the highest pressure to which oxygen 
molecules can survive, and whether oxygen forms poly- 
meric materials as found in N2 P^HTT] , CO [HI US] , and 
CO2 [iniin], and predicted in H2 [22123] • Materials un- 
der terapascal pressures are of great interest in planetary 
science, for example, the pressure at the center of Jupiter 
is estimated to be about 7 TPa [53]. Recent progress in 
dynamical shock wave [24H26| and ramped compression 
experiments [271 HH] has demonstrated that the terapas- 
cal pressure regime is becoming much more accessible. 

The use of DFT computations combined with search- 
ing methods has provided a new route for predicting 
the structures and energetics of high-pressure phases 
[131 Uni ISO]- A very recent study of phase transitions, 
melting, and chemical reactivity in CO2 found it to dis- 
sociate into carbon and oxygen above 33 GPa and 1720 
K [31], which adds further motivation for studying pure 
oxygen at high pressures. In this work we focus on oxy- 
gen at ultra-high pressures up to the terapascal regime, 
finding very surprising behavior. 

We have used the ab initio random structure search- 
ing (AIRSS) method [16, 22l [30l [32] and DFT calcula- 
tions to identify low-enthalpy structures of oxygen in the 
multi TPa range. We used the Perdew-Burke-Ernzerhof 
(PEE) [33^ generalized gradient approximation (GGA) 
exchange-correlation density functional. The searches 
were performed using the CASTEP plane-wave DFT code 
[34] and ultrasoft pseudopotentials. Searches were per- 
formed at selected combinations of 0.5, 1, 2, 3, 4, 5 and 
8 TPa, with 6, 8, 9, 10, 12 and 16 oxygen atoms per 
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cell. A total of about 14000 relaxed structures were gen- 
erated in the searches. The enthalpy-pressure relations 
were then recalculated using very hard projector aug- 
mented wave (PAW) pseudopotentials and the VASP 
code [35] with a plane-wave basis set cutoff energy of 900 
eV. Phonon and electron-phonon coupling calculations 
were performed using DFT perturbation theory and the 
ABINIT code [36], norm-conserving pseudopotentials, the 
PBE functional, and an energy cutoff of 1632 eV. Further 
details of the calculations are given in the supplementary 
material [57] . 
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FIG. 1. Enthalpy-pressure relations for solid oxygen. The up- 
per right inset shows the enthalpies of the molecular phases, 
while the insets in the lower left corner show the square spiral 
structure of Mi/acd and zig-zag chains of Cmcm, respec- 
tively. The green bubbles represent the electron lone pairs. 

The enthalpy-pressure relations for the most interest- 
ing structures are plotted in Fig. [l] We predict that at 
about 1.9 TPa molecular oxygen should transform into 
a square spiral-like structure belonging to space group 
lAi/acd, as shown in the inset in the lower left corner of 
Fig. [T] Similar structures have been found in other heav- 
ier group VI elements, such as sulfur and selenium |38| . 
The shortest 0-0 bond length in lii/acd at 2 TPa is 
about 1.15 A, which is longer than that of the molecular 
phase at the same pressure (about 1.03 A). The shortest 
0-0 distance between the chains is about 1.55 A, and the 
distance between the axes of neighboring square spirals 
is half of the lattice vector. The 0-0-0 angle in the helix 
is about 98.7°. 

At about 3.0 TPa, the square-spiral phase transforms 
into a zig-zag chain-like phase of Cmcm symmetry, as 
shown in Fig. [Tj The 0-0 bond length within the chains 
is 1.10 A at 3.5 TPa, while the inter-chain 0-0 separation 



in the plane is 1.33 A. The 0-0-0 angle within the chain 
is 102.8°. With increasing pressure the chains in each 
layer of Cmcm gradually merge and at about 9.5 TPa a 
structure of Fmmm symmetry forms in which each atom 
has four nearest neighbors at a separation of 1.05 A. 

The low-pressure molecular phases are very close in en- 
thalpy, as the inset of Fig. [T] shows. Among these phases, 
PQ^/mmc-A, C2/m-2, C2/m-8 and R3m become stable 
in turn on increasing the pressure from 0.5 to 1.9 TPa. 
More precise methods than DFT might be required to 
clarify the structural sequence as the enthalpy differences 
are so small. Although the 0=0 bond in the O2 molecule 
is not as strong as the N=N and C=0 triple bonds, O2 
polymerizes at a much higher pressure than N2 (~110 
GPa [TS]) and CO (experiments show that CO polymer- 
izes at about 5 GPa [TF while calculations predict that it 
could even polymerize at ambient pressure [19] ) . The fact 
that the molecular phases persist to pressures as high as 
1.9 TPa is intriguing. Electron counting arguments indi- 
cate that an oxygen atom can form two covalent bonds, 
either a double bond as in the O2 molecule or two single 
bonds, as would be expected in polymerized oxygen. The 
double bonds of the molecule are shorter and stronger 
than the single bonds of polymerized forms, and the lone 
electron pairs on the oxygen atoms result in bent bonds, 
as in the ozone O3 molecule, and give rise to strongly re- 
pulsive interactions between the molecules and between 
the polymeric chains. The great reluctance of oxygen to 
form more than two covalent bonds, the requirement that 
these bonds be bent, and the lone-pair repulsion greatly 
limits the forms of dense structures which may occur. 
The square spirals of the /4i /acd structure are not well 
packed and, while the zig-zag chains of Cmcm pack more 
efhciently, there is no bonding between chains so that 
very dense structures are not formed. These factors lead 
to the persistence of molecular forms to high pressures, 
and the absence of stable framework structures up to at 
least 10 TPa. 

We studied other possible candidate structural types 
such as buckled octagons of the type found in the S-I 
phase of sulfur [Fddd] and the 3- fold spiral- like P32 21 
structure of S-II [38 , but these were found to be much 
less stable than the best structures obtained in our 
searches. 

The electronic band structures and projected elec- 
tronic densities of states (PDoS) of Rim at 1.8 TPa, 
I All acd at 2.0 TPa and Cmcm at 3.5 TPa are shown 
in Fig. [2] decomposed into s and p components. As 
mentioned above, the enthalpies of the molecular phases 
are very close, but the i?3m structure is the most likely 
molecular ground state prior to polymerization at 1.9 
TPa. Our DFT calculations suggest that the R3m struc- 
ture is metallic at 1.8 TPa and the electronic density of 
states at the Fermi level mostly derives from the p elec- 
trons. 

Surprisingly, it turns out that the square spiral-like 
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FIG. 3. Electronic densities of polymeric oxygen, (a) and (b) 
74i /acd at 2.0 TPa along the [100] and [001] directions, (c) 
and (d) Cmcm at 3.5 TPa along the [100] and [001] directions. 



lAi/acd phase has a band gap of 3.0 eV at 2.0 TPa, as 
shown in Fig. |2jb). Considering that standard density 
functionals such as PBE are well-known for underesti- 
mating band gaps, lAi/acd is expected to be a wide-gap 
semiconductor. As shown in the electron density plots of 
Fig.[3]Ja) and (b), and also in the charge density difference 
in the lower left inset of Fig. [TJ the electrons are strongly 
localized. Strong covalent bonds are formed along the 
spiral chains and the projections of the electron density 
on the [100] and [001] planes show that the spirals are 
separated and not bonded to each other. 

The calculated PDoS and band gap of the Cmcm phase 
(Fig.[2];c)) show that oxygen is metallic at 3.5 TPa. The 



density of states at the Fermi level derives approximately 
equally from the s and p electrons and is smaller than 
that of the B2m structure (Fig. [2f^a)). This indicates 
that the Cmcm phase is a weak metal and supports the 
conclusion that the lone pair electrons are not very ex- 
tended. The electron density plots of Figs. [Sf^c) and (d) 
show that the zig-zag chains are not bonded to one an- 
other and that the layers are well separated. 
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FIG. 4. Phonon dispersion relations of i?3m at 1.8 TPa (up- 
per), 14,1 /acd at 2.0 TPa (lower left) and Cmcm at 3.5 TPa 
(lower right). The blue circles with magenta filling show the 
phonon linewidth at different wave vectors. 

We investigated the dynamical stability of the R3m, 
1 4,1 /acd and Cmcm structures by calculating their 
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phonon dispersion relations, as shown in Fig.|4] All of the 
structures are mechanically stable within certain pressure 
ranges. The phonons of both the molecular R3m and 
C2/m phases (shown in supplementary [37j ) are stable at 
1.0 TPa, and i?3m is also stable at 1.8 TPa. These struc- 
tures have a group-subgroup relation which may indicate 
that oxygen molecules symmetrize upon increasing pres- 
sure near the transformation to polymeric phases. The 
molecular C oxygen is a superconductor at temperatures 
below 0.6 K [IJ. We have therefore calculated the phonon 
linewidth, Eliashberg function a^Fiuj), electron-phonon 
coupling (EPC) constant A, and the logarithmic aver- 
age phonon frequency ujiog, to investigate whether the 
molecular i?3m and polymeric Cmcm phases of oxygen 
are superconducting. 

As shown in Fig. |4] (upper), the acoustic phonon modes 
and the intramolecular LO modes of RSm have rela- 
tively large linewidths close to the F point, while the 
TO modes have large linewidths along L-F. We have 
used the Allen-Dynes modification of the McMillan equa- 
tion for the superconducting transition temperature [SS] j 
Tc = w;og/1.2exp[-1.04(l-|-A)/A-/i*(l-t-0.62A)], where 
fi* is the Coulomb pseudopotential, with typical values 
of fi* between 0.10 and 0.13. The calculated isotropic 
EPC constant for the Rim structure at 1.8 TPa is about 
0.34 and varies from 2.1 to 0.6 K, which is close to the 
value observed for the C phase of about 0.6 K [T]. 

The phonon band structure of lAi/acd at 2.0 TPa is 
shown in Fig. [4] (lower left). The phonon dispersion along 
X-P and P-N is small, indicating that the bonding be- 
tween the square spirals is weak. The phonon dispersion 
relation of Cmcm at 3.5 TPa (Fig. |4] lower right) shows 
the structure to be stable and highly anisotropic. Chain- 
like and layered structures also appear in other molecular 
systems, such as CO 1121111] and Na dT]. 

We have investigated oxygen under terapascal pres- 
sures using the AIRSS technique and ab initio DFT cal- 
culations. Oxygen molecules persist up to about 1.9 TPa, 
which is almost one order of magnitude larger than pre- 
viously reported [lOj and also much larger than other 
diatomic molecules such as H2, N2 and CO. Oxygen poly- 
merizes at about 1.9 TPa giving a square spiral-like struc- 
ture of space group Mi/acd, which is a semiconductor 
with a DFT band gap of about 3.0 eV. The R3m struc- 
ture is superconducting at 1.8 TPa and has a similar Tc 
to the C phase. R3m is slightly more stable than other 
molecular forms at high pressures and its metallic oxygen 
molecules symmetrize before polymerization. At about 
3.0 TPa, the system collapses to a denser metallic Cmcm 
phase, which consists of zig-zag chains. The stability 
of molecular and chain-like phases of highly-compressed 
solid oxygen demonstrates that the electron counting rule 
is effective up to multi-TPa pressures. Merging of the 
chains in each layer of Cmcm leads to the formation of 
a phase of Fmmm symmetry at 9.5 TPa in which each 
atom has four nearest neighbors, which signals the even- 



tual failure of the electron counting rule to predict the 
bonding. The strongly repulsive interactions between the 
lone pairs hinders the polymerization of the molecules 
and the formation of more highly coordinated structures. 
The electronic structure of oxygen changes dramatically 
with increasing pressure. It transforms from the insulat- 
ing e phase to the molecular metal C phase, which persists 
up to 1.9 TPa, where it transforms into the semicon- 
ducting lAi/acd polymeric phase, eventually returning 
to metallic behaviour with the appearance of the Cmcm 
phase at about 3.0 TPa. 
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